Type-1.5 superconductivity in multiband systems: the effects of interband couplings 
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In contrast to single-component superconductors, which are described at the level of Ginzburg- 
Landau theory by a single parameter k and are divided in type-I k < l/v2 and type-II n > l/v2 
classes, two-component systems in general possess three fundamental length scales and have been 
shown to possess a separate "type-1.5" superconducting state^l In that state, as a consequence of 
the extra fundamental length scale, vortices attract one another at long range but repel at shorter 
ranges, and therefore should form clusters in low magnetic fields. In this work we investigate the 
appearance of type-1.5 superconductivity and the interpretation of the fundamental length scales 
in the case of two active bands with substantial interband couplings such as intrinsic Josephson 
coupling, mixed gradient coupling and density-density interactions. We show that in the presence 
of substantial intercomponent interactions of the above types the system supports type-1.5 super- 
conductivity with fundamental length scales being associated with the mass of the gauge field and 
two masses of normal modes represented by mixed combinations of the density fields. 



I. INTRODUCTION 

According to Ginzburg-Landau theory, a conventional 
superconductor near T c is described by a single complex 
order parameter field. The physics of these systems is 
governed by two fundamental length scales, the magnetic 
field penetration depth A and the coherence length £, and 
the ratio k of these determines the response to an external 
field, sorting them into two categories as follows; type-I 
when k < l/-\/2 and type-II when k > l/y/ffl. 

Typc-I superconductors expel weak magnetic fields, 
while strong fields give rise to formation of macroscopic 
normal domains with magnetic flusP. The response of 
type-II superconductors is completely different; below 
some critical value H c i, the field is expelled. Above this 
value a superconductor forms a lattice or a liquid of vor- 
tices that have a supercurrent circulating around a nor- 
mal core and carry magnetic flux through the system. 
Finally, at a higher second critical value, H c2 supercon- 
ductivity is destroyed. 

These different responses are usually viewed as conse- 
quences of the vortex interaction in these systems, the 
energy cost of a boundary between superconducting and 
normal states and the thermodynamic stability of vor- 
tex excitations. In a type-II superconductor the energy 
cost of a boundary between the normal and the supercon- 
ducting state is negative, while the interaction between 
vortices is repulsive^. This leads to a formation of stable 
vortex lattices and liquids. In type-I superconductors the 
situation is the opposite; the vortex interaction is attrac- 
tive (thus making them unstable against collapse into one 
large vortex) , while the boundary energy between normal 
and superconducting states is positive. From a thermo- 
dynamic point of view the principal difference between 
type-I and type-II states is the following: (i) In type- 
II superconductors the external magnetic field strength 
required to make formation of vortex excitations energet- 
ically preferred, H c i, is smaller than the thermodynam- 
ical magnetic field H ct (the field whose energy density 



is equal to the condensation energy of a superconduc- 
tor, i.e. the field at which the uniform superconducting 
state becomes thermodynamically unstable) ; (ii) In type- 
I superconductors the field strength required to create 
a vortex excitation is larger than the thermodynamical 
critical magnetic field i.e. vortices cannot form. One can 
distinguish also a special "zero measure" boundary case 
where k has a critical value exactly at the type-I/type-II 
boundary, which in the most common GL model parame- 
terization corresponds to k = 1/ y/2. In that case vortices 
do not interact^ in the Ginzburg-Landau theory. 

The above circumstances result in a situation where, in 
a strong external magnetic field, type-I superconductors 
usually have a tendency to minimize boundary energy 
between the normal and superconducting states, leading 
to a formation of large inclusions of normal phase which 
frequently have laminar structure^. 

Recently there has been increased interest in supercon- 
ductors with several superconducting components. The 
main situations where multiple superconducting compo- 
nents arise are (i) multiband superconductors^!], (ii) 
mixtures of independently conserved condensates such as 
the projected superconductivity in metallic hydrogen and 
hydrogen rich alloys^tliJ a nd (iii) superconductors with 
other than s-wave pairing symmetries. In this work we 
focus on the cases (i) and (ii). The principal difference 
between the cases (i) and (ii) is the absence of the inter- 
component Josephson coupling in case (ii). 

In two-band superconductors (i) the superconducting 
components originate from electronic Cooper pairing in 
different bands^. Therefore these condensates could not a 
priori be expected to be independently conserved. This, 
at the level of effective models should manifest itself in 
a rather generic presence of intercomponent Josephson 
coupling. 

In the case (ii) two superconducting components 
were predicted to originate from electronic and pro- 
tonic Cooper pairing in metallic hydrogen or hydrogen- 
rich alloys. In the projected liquid metallic deuterium 
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or deuterium-rich alloys, electronic superconductivity 
was predicted to coexist at ultra high pressures with 
deuteronic condensatiorP^HS] Because electrons cannot 
be converted to protons or deuterons the condensates 
are independently conserved, and therefore in the effec- 
tive model intercomponent Josephson coupling is forbid- 
den on symmetry grounds. These states are currently 
a subject of a renewed experimental pursuit. They are 
expected to arise at high but experimentally accessible 
pressures (« AOOGPa). Current static compression ex- 
periments achieve pressures of f=a 350GPa with pressures 
of an order of \TPa being anticipated in diamond anvil 
cell experiments due to the recent availability of ultra 
hard diamonds. Similar two-charged component models 
were discussed in the context of the physics of neutron 
stars where they represent coexistent protonic and E~- 
hyperon Cooper pairs in the neutron star interior^. 

This wide variety of systems raises the need to un- 
derstand and classify the possible magnetic responses of 
multicomponent superconductors. It was discussed re- 
cently that in multicomponent systems the magnetic re- 
sponse is much more complex than in ordinary systems, 
and that the type-I/type-II dichotomy is not sufficient 
for classification. Rather, in a wide range of parameters, 
as a consequence of the existence of three fundamental 
length scales, there is a separate superconducting regime 
where vortices have long-range attractive, short-range re- 
pulsive interaction and form vortex clust ers i mmersed in 
domains of two-c ompo nent Meissner state^l Recent ex- 
perimental workii^li^ have put forward the suggestion 
that this state is realized in the two-band material MgB2, 
which sparked growing interest in this topic. In particu- 
lar questions were raised over whether this "type-1.5" su- 
perconducting regime (as it was termed by Moshchalkov 
et aP^, for recent works see^) is possible even in principle 
in the case of various non vanishing couplings (e.g. in- 
trinsic Josephson coupling, mixed gradient couplings etc) 
between superconducting components in different bands. 

In this work we report a study of the appearance 
of type-1.5 superconductivity especially focusing on the 
case of multiband superconductivity, demonstrating the 
persistence of this type of superconductivity in the pres- 
ence of various kinds of intercomponent couplings (such 
as interband Josephson coupling, mixed gradient cou- 
pling, density-density, and other kinds coupling). 



A. Type-1.5 superconductivity 

The possibility of a new type of superconductivity, 
distinct from the type-I and type-II in multicomponent 
systemJ^ is based on the following considerations. In 
principle the boundary problem in the Ginzburg-Landau 
type of equations in the presence of phase winding is 
not, from a rigorous point of view, reducible to a one- 
dimensional problem in general. Furthermore, as dis- 
cussed irpJ^I, in general in two-component models there 
are three fundamental length scales which renders the 
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FIG. 1: A comparison of the magnetic phase diagrams of 
clean bulk type-I,type-II and type-1.5 superconductors at 
zero temperature. The semi-Meissner state is a macroscopic 
phase separation into two-component Meissner state and vor- 
tex clusters where one of the density modes is suppressed by 
core overlaps. 



model impossible to parametrize in terms of a single di- 
mensionless parameter k. In the case where the conden- 
sates are not coupled by interband Josephson coupling 
but only by the vector potential these length scales are 
the two independent coherence lengths (set by the in- 
verse masses of the corresponding scalar density fields) 
and magnetic field penetration length (set by the in- 
verse mass acquired by the gauge field). In contrast, 
in the case where the condensates are coupled by inter- 
band Josephson terms, one cannot distinguish indepen- 
dent coherence lengths attributed to different conden- 
sates. Nonetheless, in this case the density variations 
can also possess two fundamental length scaleiP, in con- 
trast to single-component theories. We elaborate on this 
fact below. In 12 vortex solutions in two-component the- 
ories were found which have non-monotonic vortex inter- 
action, with a long range attractive part determined by 
a dominant density-density interaction and a short range 
repulsive part produced by current-current and electro- 
magnetic interactions. An important circumstance which 
was demonstrated was that these vortices are thermody- 
namically stable in spite of the existence of the attractive 
tail in the interaction. 

A non-monotonic intervortex interaction potential 
should result in the formation of vortex clusters in low 
magnetic field immersed into the vortexless areas, a state 
referred to irP as the "semi-Meissner state" . Figure [l] 
shows the schematic phase diagram of a type-1.5 super- 
conductor. 

If the vortices form clusters one cannot use the 
usual one-dimensional argument concerning the energy 
of superconductor-to-normal state boundary to classify 
the magnetic response of the system. First of all, the 
energy per vortex in such a case depends on whether a 
vortex is placed in a cluster or not: i.e. formation of a 
single isolated vortex might be energetically unfavorable, 
while formation of vortex clusters is favorable, because 
in a cluster where vortices are placed in a minimum of 
the interaction potential, the energy per flux quantum is 
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smaller than that for an isolated vortex (thermodynam- 
ically the nonmonotonic two- vortex interaction potential 
predicts that the smallest energy per flux quantum will 
be in the case of a uniform lattice with spacing equal to 
the minimum of two-body intervortex potential). 

Thus, besides the energy of a vortex in a cluster, there 
appears an additional energy characteristic associated 
with the boundary of a cluster. In other words, in this 
situation, to determine the magnetic response of a system 
it is not sufficient to study the one-dimensional bound- 
ary problem nor the single-vortex problem, in contrast 
to single component systems. Moreover, in a cluster 
the system tends to minimize the boundary energy of 
a cluster (similarly to type-I behavior), while breaking 
into a lattice of one-quantum vortices inside the cluster 
(similarly to type-II systems with negative interface en- 
ergy). Thus, in an increased magnetic field the vortices 
form via a first order phase transition. A magnetic phase 
distinct from the vortex and Meissner states which then 
arises is a macroscopic phase separation into domains of 
two-component Meissner state and vortex clusters where 
one of the density modes is suppressed by core overlap. 
We summarize the basic properties of type-I, type-II and 
type-1.5 regimes in the table |TJ 

The existence of thermodynamically stable type-1.5 su- 
perconducting regimes ultimately depends on the exis- 
tence of a nonmonotonic intervortex interaction poten- 
tial. It is an important question how generic this effect is. 
In this work we mainly focus on multiband realizations 
of multicomponent superconductivity and investigate the 
effects of interband Josephson coupling, mixed gradient 
coupling, and density-density coupling terms on vortex 
interactions in two band superconductors. We show that 
(i) when these couplings are present, the system still can 
possess three fundamental length scales, in contrast to 
the two length scales in the usual single-component GL 
theory; (ii) non-monotonic interaction is possible in a 
wide parameter range in these models. 

The structure of this paper is as follows: In section II 
we introduce the model. 

In section III we present a linear theory of asymptotics 
of the vortex fields in a superconductor with two bands 
with various interband couplings. 

We begin section III by demonstrating that for a gen- 
eral form of the effective potential in a two-band (or more 
generally two-gap) Ginzburg-Landau free energy, the lin- 
ear theory gives, under quite general conditions, two fun- 
damental length scales of the variations of the densities. 
From the linearized theory we calculate the long-range in- 
tervortex interaction potentials using the two-component 
generalization of the point-vortex method^ and show 
how the non-monotonic intervortex interaction potential 
arise from the interplay of two fundamental length scales 
of the superfluid density variations and the magnetic field 
penetration length. The central point of this part is how 
the fundamental length scales are defined in the presence 
of interband coupling as well as the occurrence of "mode 
mixing" . Next we move to quantitatively study the ef- 



fects of several kinds of intercomponent couplings which 
quite generically arise in two-component theories. 

In section III (d) we demonstrate that that mixed gra- 
dient coupling can lead under certain conditions to an 
increase in the disparity of the characteristic scales of 
the density variations. 

In Section IV we present a large scale numerical study 
of the full nonlinear problem of the interaction between 
a pair of vortices. 

II. THE MODEL 

A. Free energy functional 

We study the type-1.5 regime using the following two- 
component Ginzburg-Landau (TCGL) free energy func- 
tional. 

F = ^ViXW + \(D^ 2 )(Di/j 2 y (1) 

-vRe{{Dip x ){D^)*) + ~(V x A) 2 + F p 

Here D = V + ieA, and ip a = |i/> a |e iea , a = 1,2, rep- 
resent two superfluid components which, in a two-band 
superconductor correspond to two superfluid densities in 
different bands. The term F p can contain in our analysis 
an arbitrary collection of non-gradient terms. 

A particular form of two-component GL model which 
was microscopically derived irPHSl f or two-band super- 
conductors is: 

F=ipVi)WO* + ^(£»^a)P^)' (2) 
-uRe{(D^)(Dij 2 )*} + ~(V x A) 2 

+ai|V>i| 2 + ^ilV>i| 4 +a2^ 2 | 2 + ^IV^I 4 
-^ilVill^l cos(#! - 9 2 ) + ^IViH^I 2 - 

The first two terms represent standard Ginzburg- 
Landau gradient terms, the second term represents mixed 
gradient interactions which were shown to originat e in 
two-band superconductors from impurity scattering^!. 
The next term is the magnetic field energy density and 
the remaining terms represent an effective potential. 
Here we note that a± and a 2 can invert sign at dif- 
ferent temperatures. The regime where ot\ is positive 
while a 2 is negative corresponds to the situation where 
one of the bands has no superconductivity of its own 
but nonetheless bears some superfluid density due to in- 
terband Josephson tunneling, which is represented here 
by the term ^iI^iHV^I cos(#i — 9 2 ). The type-1.5 be- 
haviour in this regime was studied irP. In this work we 
will mainly focus on the situation where both bands are 
active, i.e. olx,2 < 0. For generality we also add a higher 
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single-component Type-I 


single-component Type-II 


multi-component Type-1.5 


Characteristic 
lengths scales 


Penetration length A & coher- 
ence length £ (| < i) 


Penetration length A & coher- 
ence length £ (| > i) 


Two characteristic density varia- 
tions length scales £i,£2 and pene- 
tration length A, the nonmonotonic 
vortex interaction occurs in these 
systems typically when £i < \/2A < 


Intervortex inter- 
action 


Attractive 


Repulsive 


Attractive at long range and repul- 
sive at short range 


Energy of su- 
perconduct- 
ing/normal state 
boundary 


Positive 


Negative 


Under quite general conditions 
negative energy of superconduc- 
tor/normal interface inside a vortex 
cluster but positive energy of the 
vortex cluster's boundary 


The magnetic field 
required to form a 
vortex 


Larger than the thermodynam- 
ical critical magnetic field 


Smaller than thermodynamical 
critical magnetic field 


In different cases either (i) smaller 
than the thermodynamical critical 
magnetic field or (ii) larger than 
critical magnetic field for single vor- 
tex but smaller than critical mag- 
netic field for a vortex cluster of a 
certain critical size 


Phases in external 
magnetic field 


(i) Meissner state at low fields; 

(ii) Macroscopically large nor- 
mal domains at larger fields. 
First order phase transition be- 
tween superconducting (Meiss- 
ner) and normal states 


(i) Meissner state at low fields, 

(ii) vortex lattices/liquids at 
larger fields. Second or- 
der phase transitions between 
Meissner and vortex states and 
between vortex and normal 
states. 


(i) Meissner state at low fields 

(ii) "Semi-Meissner state": vor- 
tex clusters coexisting with Meiss- 
ner domains at intermediate fields 

(iii) Vortex lattices/liquids at larger 
fields. Vortices form via a first or- 
der phase transition. The transition 
from vortex states to normal state 
is second order. 


1 ' 1 ' / TV T \ (' 

Energy H^INj oi 
N-quantum ax- 
ially symmetric 
vortex solutions 


N < j v _ 1 tor all IN. Vor- 
tices coalesce onto a single N- 
quantum megavortex 


N > iv-1 tOT a11 N - 

N-quantum vortex decays into 
N infinitely separated single- 
quantum vortices 


There is a characteristic number N c 
such that 3£2 < E <£-» for N < 

N c , while E ^ > for N > 
N c . N-quantum vortices decay into 
vortex clusters. 



TABLE I: Basic characteristics of bulk clean superconductors in type-I, type-II and type-1.5 regimes. Here the most common 
units are used in which the value of the GL parameter which separates type-I and type-II regimes in a single-component theory 
is K a = l/y/2. 



order density-density coupling term r\ 2 \ i>\ 1 2 1 | 2 ■ We also 
consider the case of independently conserved condensates 
where the third and ninth terms in ^ are forbidden 
on symmetry grounds, that is, v = rji = (see also 
remarlPil). The equivalence mapping between our units 
and the standard textbook units is given in Appendix [A"| 
A microscopic derivation of the TCGL model (J2j) re- 
quires the fields \ip a \ to be small. However it does not in 
principle require a a to change sign at the same temper- 
ature. Moreover, as in the case of single-component GL 
theory, we expect the model ^ to give in many cases 
a qualitatively acceptable picture in lower temperature 
regimes as well. In fact, our analysis can in some cases 
give a qualitative picture for the case where one of the 
fields does not possess a GL-type effective potential be- 
cause the regime where one of the bands is in a London 
limit (i.e. it does not possess a GL effective potential but 
has a small vortex core modeled by a sharp cutoff) can be 



recovered from our analysis as a limiting case. As will be 
clear from the analysis below that regime also supports 
type-1.5 superconductivity. 



B. Basic properties of the vortex excitations. 

The only vortex solutions of the model ^ which have 
finite energy per unit length are the integer A^-flux quan- 
tum vortices which have the following phase windings 
along a contour I around the vortex core: f, V#i = 
2irN, §. V6* 2 = 2ttN. Vortices with differing phase wind- 
ings carry a fractional multiple of the magnetic flux quan- 
tum and have energy divergent with the system size. 
These solutions were investigated in detail irp3. 

In what follows we investigate only integer flux vortex 
solutions which are the energetically cheapest objects to 
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produce by means of an external field in a bulk super- 
conductor. 



III. VORTEX ASYMPTOTICS 



equation system 



(T) 



(U1,M 2 ,0) 



1 



a , -e(l + ea)(/ 1 2 + /|) = 0. 



(8) 



The key to understanding the interaction of well sepa- 
rated vortices is to analyze the large r asymptotics of the 
vortex solution. We will analyze this problem in the con- 
text of a general TCGL model whose free energy takes 
the form 



The solution we require, the vortex, has boundary be- 
haviour f a (r) — > u a , a(r) —> — 1/e as r — > oo. So, for 
large r, the quantities 



1 



e (r) = f a (r) - u a , a(r) = a(r) + - 



(9) 



1 



(5^2 - d 2 A 1 f + F p 



(3) 



where F p contains all the non-gradient terms (in particu- 
lar, but not restricted to, Josephson and density-density 
interaction terms). This free energy is consistent with 
([2]) in the case v = 0. We will show in section HID how 
to handle mixed gradient terms. The precise form of F p 
is not crucial for our analysis in this section. By gauge 
invariance, it can depend on the condensates only via 
IV'iL I ^2 1 and (if the condensates are not independently 
conserved) on 61 — 62- We will assume that F p takes 
its minimum value (which we normalize to be 0) when 
|V>i| = ui > 0, \i/j 2 \ = u 2 > and 61-62 = 0. So, either 
there is no phase coupling (F p is independent of 61 — 62) 
and the choice of 6\ — 62 = is arbitrary, or the phase 
coupling is such as to encourage phase locking. (Note 
that the case of phase anti-locked fields can trivially be 
recovered from our analysis by mapping ip2 <-> —ipi)- 

The field equations are obtained from F by demanding 
that the total free energy E — J Fdxidx2 is stationary 
with respect to all variations of ipi , ip2 and Ai . A routine 
calculation yields 



BF 

dtp* 



(4) 



(5) 



This triple of coupled nonlinear partial differential equa- 
tions supports solutions of the form 



(Ai,A 2 ) 



fa(ry e 
a(r) 



(— sinfl, cos( 



(6) 



where /i,/2,a are real profile functions. Note that in 
some cases mixed gradient terms favour non-axially sym- 
metric solutions. In this section we consider only axially 
symmetric vortices. Fields within the above ansatz sat- 
isfy the field equations if and only if the profile functions 
fi( r )i f2(r), a(r) satisfy the coupled ordinary differential 



are small and so should, to leading order, satisfy the lin- 
earization of ([7]),(|8| about (ui, u%, — 1/e). That is, at 
large r, 



1 , 



^ I 2/ 2 1 2\ 

-a — e (u 1 + u 2 )a 



b=l 





(10) 



(11) 



where T-L is the Hessian matrix of F p (|^i|, l^l) 0) about 
its minimum 



T~Lab 



d 2 F n 



d\lPa\d\TP b 



(12) 



(U1,M 2 ,0) 



So a decouples from ei , £2 asymptotically, and we see 
immediately that 



a(r) = qorK^Ar), fi A 



e\j u\ + u 



(13) 



where K n denotes the nth modified Bessel's function of 
the second kincP3, and qo is an unknown real constant. 
Hence, at large r, 



— + q^Ki {fiAr) ) (— sin 6, cos ( 

er 



Since, for all n, 

K n (s) 



00, 



(14) 



(15) 



it follows that the magnetic field decays exponentially as 
a function of r, with length scale (penetration depth) 



A : 



1 



1 



V M ? + ■ 



(16) 



By contrast, (10) represents, in general, a coupled 



pair of ordinary differential equations for £1,62- Since 
(til, 1/2,0) is a minimum of -f n p (|'0i | , |'02 1 1 ^1 — $2), the 
Hessian matrix % is a positive definite symmetric 2x2 
real matrix. Hence its eigenvalues, [if, sa Y: are rea l 
and positive, and its eigenvectors, v±, v 2 say, form an or- 
thonormal basis for K 2 . Expanding e = (ei,£2) T in the 
basis t?i, t»2 



e M = Xi( r ) v i + X2(r)v 2 , 



(17) 
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we see that \i > Xi satisfy the uncoupled pair of ordinary 
differential equations 

X'a + \x' a = tiXa, (18) 

whence 

Xa(r) = q a K Q {ji a r) (19) 

for some (unknown) constants q\ , q 2 . Since V\ , v 2 are or- 
thonormal, there is an angle O, which we call the mixing 
angle, such that the eigenvectors of H are 

Hence, at large r the density helds behave as 

ipl ~ [ui + qi cos QKo^r) - q 2 sin QK (^i 2 r)]e lB 
ip 2 ~ [u 2 + q 1 smQK (p, 1 r) +q 2 cosQK (fi 2 r)]e i9 21) 

where, once again, Kq is a Bessel function. 
From this analysis it follows that: 

1. In general there are three fundamental length scales 
in the problem (in contrast to the two length 
scales of one-component Ginzburg-Landau theory) 
which manifest themselves in the vortex asymp- 
totics, namely and 1/^2- 

2. These are constructed from the vacuum expecta- 
tion values u a of \ip a \ (in the case of I /[J. a) and 
from the eigenvalues of H, the Hessian matrix of 
F p about the vacuum (i.e. the ground state). 

3. 1/fJ-A can be interpreted as the London penetration 
length of the magnetic field. 

4. However, unless the mixing angle is a multiple 
of 7r/2, 1/fii and 1/^*2 cannot be interpreted as 
the coherence lengths of if>i,ip2 in the usual sense. 
This is because the normal modes of the field theory 
close to the vacuum are not \tp a \ — u a , but rather 

Xi = {\ipi\ -ui)cos8- (|V>2| - U2)sinG 
X2 = Qipi\ -ui)sin6- (\tp 2 \ -U2)cos6 

obtained by rotating through the mixing angle O, 
which is also determined by %. Therefore in general 
(e.g. in the presence of intercomponent Josephson 
coupling) for a one-flux quantum axially symmetric 
vortex, the recovery of both fields ip a at very long 
range will be according to the same exponential 
law, set by the smaller of the masses Hi,n 2 ; One 
should use the representation in terms of the fields 
Xi 2 to be handle properly the two length scales 
associated with the density recovery. 



5. This analysis tells us only about the vortex struc- 
ture at large r. It gives no direct information on 
the vortex core, which is important to understand 
quantitatively the nature of the vortex interactions 
at intermediate and short distances which will be 
studied numerically in the section V. 

Since the gauge field mediates a repulsive force between 
vortices, while the condensate fields mediate an attrac- 
tive force, it is clear that we can read off from the above 
analysis the condition under which the intervortex force 
is attractive at long range: we require that 1 / fiA is not 
the longest of the three length scales, or, more explic- 
itly, that (at least) one of the eigenvalues of H should be 
less that fi\ = e 2 (u\ + u 2 ). We can predict an explicit 
formula for the long range two-vortex interaction poten- 
tial, using the point vortex formalism^ ( a brief review of 
the method in given in Appendix [C| . This rests on the 
observation that, far from its core, the fields of the vor- 
tex are identical to those of a hypothetical point particle 
in a linear theory with two Klein-Gordon fields (xi and 
Xi above) of mass and a vector field (A) of mass 

fiA- The point particle carries scalar monopole charges 
2irqi and 2irq 2 and a magnetic dipole moment 27rgo- Two 
such hypothetical particles held distance r apart would 
experience an interaction potential 

V(r) = 2tt [qlK (fi A r) - fiR^r) - q 2 2 K (n 2 r)} . 

(22) 

This formula reproduces the prediction explained above: 
the long range interaction will be attractive if (at least) 
one of fii, iu>2 is less than /i^. 

One can ask, retrospectively, whether the ap- 
proximation of linearizing in the small quantities 
a { r )i Xi( r )> X2( r ) is well justified. Rigorous analysis of 
the single component modeP^l shows that if either of the 
scalar mode masses, fj, 2 say, exceeds 2/x^, then quadratic 
terms in a become comparable at large r with linear 
terms in X2, so that the equation for X2 should include 
extra terms. In this case, X2 decays like K {^Af) 2 rather 
than Ko(/j, 2 r). One should note, however that, unless 
/ii > 2/a a also, the leading term in ( |21[ ), decaying like 
K (pir), is still correct, and it is only the leading term 
which determines the nature (attractive or repulsive) of 
the intervortex interactions at long range. The case of 
interest to us is when the long-range force is attractive, 
that is, when at least one of \i\,ii 2 is less than (1a, so 
the linearized analysis presented above suffices for our 
purposes. 

A. The U(l) x 17(1) symmetric model 

We first illustrate the above analysis in the case of 
U(l) x U(l) condensates coupled only by a gauge fielcP 
where 

F p = a!|Vi| 2 + y |V>i| 4 + aM 2 + y> 2 1 4 + constant, 

(23) 
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with both at < and ct2 < 0. Here F p is independent of 
9 1 — 62 and its minimum occurs at \ip a \ = u a where 



Pa 



The Hessian matrix of F p at (u±, 112, 0) is 



V. = 




(24) 



(25) 



whose eigenvalues are fj, a = — 4a a , with corresponding 
eigenvectors v\ = (1,0) , v% = (0, 1) T . Hence, the mix- 
ing angle is = 0, the penetration depth is 



A = — = e 

(J'A 



-a 2 



Pi P2 



(26) 



and the density decay lengths are the usual coherence 
lengths 



1 



fj, a 2^-a a 



(27) 



The criterion for long-range vortex attraction therefore 
amounts to the requirement that one of the coherence 
lengths is larger than the magnetic field penetration 
length, 



iin{2- v / — ai, 2y/— a 2 } < e 



-Ol 

Pi 



~«2 
P2 ' 



(28) 



Note this criterion indicates only a long range attrac- 
tion. For a realization of the type-1.5 regime one should 
additionally require short range repulsion and thermody- 
namic stability of a vortex cluster (i.e. a vortex cluster 
should become energetically favorable to form in exter- 
nal fields smaller than the thermodynamical critical mag- 
netic field). 

Note also that, in systems with independently con- 
served condensates, such GL models arise from expansion 
of the free energy in powers (1 — T/T cl ) and (1 — T/T c2 ) 
near two critical temperatures (such systems are, for ex- 
ample, the projected superconducting states of metallic 
hydrogen or condensate mixtures in neutron stars). In 
general T c i and T C 2 can be quite different making two 
such expansions at the same temperature formally im- 
possible. In that case the more suitable model is a Lon- 
don model for one component (i.e. |t/>i| const except a 
core cutoff) coupled to a GL model of the second compo- 
nent. When it is the component with "short" coherence 
length which is modeled by the London limit, we recover 
a description of the type-1.5 behavior in that case from 
our above analysis as a simple limit. 



B. Josephson coupling 

We next consider how this picture is influenced by the 
addition of an interband Josephson term which breaks 



the U(l) x U(l) symmetry to U(l) (note that this in 
particular implies that there is only one critical temper- 
ature), 



F p - 771^1 1 1^2 1 cos(0i 



(29) 



where F p is the free energy defined in (23) and the 



Josephson coupling 771 > 0, so that F p is minimized when 
B\ — 62 = 0. Adding this term changes the vacuum ex- 
pectation values u a of the fields. To find 1*1,1*2 we must 
solve 



dF n 



dF r , 



that is, 



1ct\Ui 
2a 2 U2 



2Piu\ 
2p 2 u 2 2 



0. 



V1U2 

m u i- 



(30) 



(31) 
(32) 



Unfortunately, it is not possible to solve these equations 
explicitly, except in special cases. For particular values 
of the parameters a. a , p a , rji they can easily be solved nu- 
merically, as can the eigenvalue problem for %. Note 
that like in the case of uncoupled bands, there are in 
general three fundamental length scales also in the pres- 
ence of the Josephson term which can then be computed. 
We present numerical analysis of this problem in the full 
Ginzburg-Landau model in section |IV| To make analyt- 
ical advance in this section we treat the 771 dependence 
of the length scales perturbatively. That is, we will con- 
struct Taylor expansions for it a (?yi), ^(771), n a (Vi) an< A 
0(?7i). To keep the presentation simple, we will work to 
order rj\, so the results will give the leading correction 
to the formulae of the previous section as the Josephson 
coupling 771 is "turned on" . Higher order corrections are 
easily computed but do not give much extra insight. 

Let us denote quantities defined in the uncoupled 
model (at 7/1 = 0) with a hat, so u a = y/—ot a /Pa are 
the uncoupled vacuum expectation values of \ip a \, for ex- 
ample. Let it(77i) = (mi(t7i), M2(t7i)) t and 



\ dF p /d\ip 2 \ 



(33) 



Then, by definition G (11(771)) = for all 771, and G(u) = 0. 
Differentiating with respect to 771 (denoted by a prime), 
we see that 





u'(0) 



dG 

dm 



(u) + Hu'{0) 



~, idG,^ 
drji 



Ar 2 








(34) 



V A2 

Hence the ground state densities receive a correction lin 



ear m r\\ 



i \ /. M 2 . 9\ 

ui{r]i) = ui + ^rj! + O^) 
Mi 

m 2 (?7i) = "2 + + 0{rjl). 

M 2 



(35) 



From this expression for ground state densities one can 
readily calculate the gauge field mass [j,a{vi)i whose in- 
verse gives the London penetration length. One sees that 



1 

IP- 
Mi 



i)»?i + o(vr). w) 

M 2 , 



So the effect of the Josephson coupling is always to in- 
crease the vacuum expectation values of \ip a \, and hence 
to decrease the penetration depth 1 / fiA ■ 

The other two length scales are the eigenvalues of T-L 
where 



H 



if 



<it, 



d\lp a \d\lP b 
Tiab + Vl 



(F p -m\fl>i\\H) 



d 3 P v 



d\^ a \d\i> b \d\i> c 
-mil -s ab ) + 0^1). 



u' c (0) 



Now 



d 3 F n 



12/3 a |^ a | ifa = 6 = c 
otherwise 



(37) 



(38) 



n = u + m 



12/3iuiu'i(0) 



12/3 2 u 2 m' 2 (0) 



0{r,l) 



fl\ + 377iM 2 /U! 



-Vi 



A| + 3?7iUi/m 2 



0{vl). (39) 



So the eigenvalues A = p? a satisfy the characteristic equa- 
tion 



A! + 3r/i ^ - A J ^ + % ^ - A J + 0(r]f) = 0. 

(40) 

Hence 



m? = tf + ^ + ofa?) 
"1 



^ = A 2 +3??1 "1+0(^2) 
"2 



(41) 



As with the penetration depth, the effect of Josephson 
coupling (at leading order), is to decrease the character- 
istic length scales l//i of both normal modes. Another 
effect of Josephson coupling is mixing the fields. Let us 
now compute the corresponding mixing angle 0(rji). Re- 
call that this is, by definition, the angle such that 



v(Vi) 



cos 6(771) 
sin Q(r/i) 



(42) 



is the eigenvector of T-L with eigenvalue /if. We know that 
v(0) = v= (1,0) T , that v(t]i) ■ v(r]i) — 1 and that 

M(rn)v(rn) = 0, where Mfa) = H(r n ) - Mvifh 

(43) 

-o{vi) 

(44) 



for all 771. As computed above, 



1 



Differentiating (43) and v{rj\) ■ v(rji) — 1 with respect to 
771 yields 

M'{0)v + M(0)u'(0) = 0, v(0) ■ v'{0) = (45) 
which can be solved for v'(0). One finds that 



v(Vi 







Vi 



(Al 





-A?) 



Hence, the mixing angle is 
®(Vi) 



Vi 



Al 



Ml 



0{ V l). (46) 



(47) 



Thus the Josephson term produces mode mixing. 
Clearly, this perturbative expansion is well-defined only 
if Ai 7^ A2- In the case where Ai = A2, Ti = AP2 and 
the assertion that v(Q) — (1,0) T is arbitrary: any or- 
thonormal pair of vectors can be taken as the eigenvec- 
tors associated with A?, Ai- Hence, the notion of "mixing 
angle" is ill-defined in this case, and is not amenable to 
perturbative calculation. 

The normal modes are associated with the following 
combinations of the \ip a \ fields 



Xl = (|V>l| - Ml) COS 

Xi = (|V>i| - Mi) sin 



Vi 



A 1 * 9 

M ~ Mi 
Vi 

A2 - A? . 



- (|^2 1 - M 2 )sin 
(l^al ~ m 2 )cos 



A2-A1. 

Vi 

A2 - A?. 



Thus one can associate "coherence lengths" of these fields 
with the inverse masses of the model ( |41[ ) which are func- 
tions of the coherence lengths (Aa 1 ) an d vacuum field 
densities (u a ) defined in the Josephson- uncoupled the- 
ory, and the strength of the Josephson coupling i]i . Note 
that, returning to the original fields \ipa\, the very long- 
range behavior of both of these density fields is governed 
by whichever of the fields xi,2 has the slower recovery 
rate. This implies that at very long range both fields 
\ip a \ should have the same exponential recovery law set 
by the smaller of /i a . The physical meaning of mode 
mixing is that the variation of the original density fields 
\ip a \ acquires two length scales and one should rotate the 
fields through the mixing angle O to determine the nor- 
mal modes xi j2 whose recovery rates are governed by 
different single exponential laws. 

Another point to note here is that, from a quantita- 
tive point of view, turning on a small Josephson coupling 
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does not radically alter the integer flux vortices. For 
small 7/1 , there is a small correction to each of the length 
scales (all three length scales become smaller), and there 
is a small amount of normal mode mixing (measured by 
0(?7i)). Therefore for composite integer flux vortices the 
addition of Josephson coupling does not represent any 
kind of singular perturbation. 



1. Comparison with the case of passive second band 

It is interesting to compare these results with the case 
where one of the bands, ip2 say, is passive, and has su- 
perconductivity only by virtue of the Josephson coupling 
terrrP. In this case, the free energy F p has ct2 > 0. 
In the uncoupled model (at rji — 0), F p is minimized 
when |?/>i| = ill = \] —u\j fix and — U2 — 0. The 
gauge field has mass {1a — eui, there is no mode mix- 
ing, and the condensates have masses fax — 2y/—ax and 
{±2 = \J1ol2- The vortex solution has ^2=0 everywhere 
and is identical to the Abrikosov vortex of the single com- 
ponent GL theory with ifj2 set to zero. As the Josephson 
coupling is turned on, ip2 acquires a vacuum density of 
order r)\, mode-mixing develops, and the three length 
scales acquire corrections. Repeating the arguments of 
the a?2 < case, one finds that 



refer to the uncoupled 772 = model): 



"2(771) 

Ma(?7i) 2 
M?7i) 2 

9(77!) 



U\ + 
Ui 

T2V1 

V-2 

Aa 
A? - 

£2 " 



+ 0(4) 

-o(vj) 
o(vl) 

0(r,l) 



Al 



- 2 
Mi 



(48) 



A significant difference from the case of two active bands 
(«2 < 0) is that all three of the length scales receive 
corrections only at order rfc. Nevertheless, there is mode- 
mixing at order 77}. 



C. Density-density coupling 

A similar perturbative analysis of the case when there 
is bi-quadratic density-density coupling 



(49) 



can be carried out. Since the calculations are similar, we 
merely record the results (once again, hatted parameters 



ux{m) = + — T2-V2 + 0{r) 2 ) 
Mi 

U2 + -^V2 + 0{T)l) 



U2(?72) 

Ma(t? 2 ) 2 



M2M 2 



A2 

{i\ + 2u\ul 



Ml 



M2 , 



- 9 

Mi 



- 9 
M2 



\% ) ufa + 0(i£) 
Mi , 

'"2 ) ^>i»i) 

M2, 



2miU 2 . n( 2n 

-To T2 ? 72 + 0(77 2 ). 

Ml M2 



(50) 



The effect of the extra term is to reduce (if 772 > 0) or in- 
crease (if 772 < 0) all three length scales and to introduce 
mode mixing. We present numerical analysis of this kind 



of coupling in section IV 



D. Mixed gradient terms 

In this section we consider the case where the free en- 
ergy has gradient-gradient coupling terms, 



F = ~(D i ip 1 )*D i ip 1 + ~(D i ij; 2 )*D i ii; 2 

--[(AWAlk + (A^)*A^i]+^ (51) 



where Fp(\tpi\,\ip2\,8i — 8 2 ) is, as before, a non-negative 
function minimized at (ui,tt2,0). In contrast to the pre- 
vious two cases, this we can treat exactly, without resort- 
ing to power series expansion in the coupling parameter 
v. We can assume v > without loss of generality (the 
case v < is obtained by mapping t/> 2 i-> —1P2), and we 
must have v < 1, or else F is not positive definite. 

This case does not fit into the general analysis pre- 
sented above. Nonetheless a similar method, with some 
modification, can be applied. 

The field equations are 



A A [i>2 - wpi) 



dr 2 



(52) 
(53) 



di(diAj - djAi) = eIm(^iA(V'i ~ ^2) 

+ VIA'(^2 - WO) (54) 

which support vortex solutions of the form (|6| provided 
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the profile functions obey the coupled system 



dr 2 



1 d 

r dr 



1 



(1 + ea) 




= I dF p /d\il>x\ 
\ dF p /d\ih\ 



a" - -a' - e(l + ea)(f 2 - 2vfJ 2 + f 2 ) = 0, 
r 



where 



P = 



(/i,/ 2 ,0) 

(55) 
(56) 

The vortex boundary conditions are f a (r) ~ > u a and 
a(r) — > — 1/e as r — > oo. Note that these are indepen- 
dent of v. We again define e(r) = (fi(r) — u\, /2(f) — u 2 ) T 
and a(r) — a(r) + 1/e, and linearize the system about 
£i = £ 2 = a = 0: 




dr 2 



1 d 

r dr 



Pe = Ue 



l -a> 
r 



e 2 (u 2 — 1vu\Ui + u\)a 



0. 



(57) 



We immediately see that A behaves asymptotically as in 
(fl4|, but with 



Ha = ey u i — 2vu\u 2 



(58) 



The effect of the gradient-gradient coupling is thus to 
increase the penetration depth 1/ha- Note that the ef- 
fect would be opposite if there is a competing sufficiently 
strong Josephson term which enforces phase anti-locking 
in the vacuum, that is phase difference 6\ — 8 2 = n. 

To decouple the pair of equations for ei,e 2 we must 
expand e in a basis of eigenvectors, not of 'H, but rather 
of 



%{y) = P(^)- 1 'H. 



(59) 



Note that this matrix is not in general symmetric. 
Nonetheless, it can be shown that its eigenvalues are real 
and positive for all < v < 1 (see the appendix B). 
Let Hi be the eigenvalues of H and v±, v 2 be the cor- 
responding eigenvectors. Then the condensate fields at 
large r take the form 

+ qiK ((iir)vi + q2l<o(p-2r)v 2 ^j e lB 

(60) 

Once again, the condition for long-range attraction is 

min{nl,(4} < n\ (61) 

but now Hi,n 2 are the eigenvalues of P(v)^ 1 H, not H, 
and ha depends on U\, u 2 and v : as in (158]). Note that 




h\h\ = detn = det^/detP(r7i) 



detn 

1 - v 2 



(62) 



so the effect of the coupling must be to increase (at least) 
one of Hi,2 and hence to decrease (at least) one of the nor- 
mal mode recovery length scales. Since H is not symmet- 
ric, there is no reason why v±, v 2 should be orthogonal, 
so it is not possible to define a single mixing angle in this 
case. 

To illustrate, consider the simplest case, where F p is 
defined as in ( 23 1 . Then u a = \f—a a /Pa and the uncou- 
pled v = model has fi 2 — —^a a and £l 2 a — e 2 {u\ + u\). 
For v > 0, the vacuum expectation values do not change, 
but the penetration depth increases, since 



h\ = P>a~ ^VUyU 2 . 



(63) 



Furthermore 



H(y) = P(yy 



i / Hi v t4 



A 2 



V 2 \ Vjj 2 jj 2 



whose eigenvalues are 



(64) 



i 



2(l-i/ 2 ) 



, , i A? + Al ± v (A! - Ai) 2 + 4^ 2 A 2 Ai 



(65) 



Now Hi 2 i. v ) are the new fundamental length scales which 



control the variation of the density fields ( 60 1 . Without 



loss of generality, we may assume that Hi > Hi (if Hi < 
Hi then we simply swap the labels of the condensates). In 
this case, for v > it is clear from the above expression 
that 



Mi > 



(T ^y (a? + As + \/(Af-A!) 2 ) 

so when < v < 1, Hi{ v ) > Al- Recall that 
/i 2 ^2 = det % 



A? 



AfAj 
i-i/ 2 



(66) 



(67) 



so 



A? 



Ai 



M f l 



< Ai 



(68) 



by (66l, and hence Hi(y) < A2 when < f < 1. Ira £/us 



case, £/ie effect of 'gradient- gradient coupling is to decrease 
the smaller of the normal mode decay lengths, Hi 1 ; an d 
increase the larger, h 2 X ■ Thus gradient coupling tends to 
increase the disparity in these length scales. 

As in sections IIIIBI and MI CI it is instructive to see 
what happens to the parameters of the uncoupled model 
(y = 0) as v > is turned on. This can be extracted 
from the above formulae by expanding in v, keeping only 
terms up to order v 1 . One sees that 



U a (v) 



u a exactly 

Ha ~ c Zvu\u 2 exactly 

A2 + <V) 



(69) 
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so, to leading order in v, the only length scale that 
changes is the penetration depth 1/fiA, which increases. 
The eigenvectors of T-L are 



1 



-0(f 2 ), v 2 = \^W^F 2 j + 0{is 2 

(70) 

which, one should note, are not orthogonal: the angle 
between them is f — v-\-0(v 2 ). It follows that the vortex 
has asymptotic densities (at large r) 

IViI ~ fii + gi^o(Air) - g 2 2A ^ 2 ifo(A20 + Ofr 2 ) 

Ml l L 2 

- £ 2 + q 2 K {(i 2 r) + ^^ Koifar) + 0(i/frl) 
Mi M2 

where 91,92 are unknown constants. So, while the "co- 
herence lengths" remain unchanged to leading order, the 
normal modes with which they are associated do receive 
a correction at order v . 

Finally, we remark that an alternative approach to 
handling gradient-gradient terms is to remove them from 
F from the outset by a linear redefinition of the fields: es- 
sentiallyone expands (ipx, V>2) T in a basis of eigenvectors 
of P{v)^i. This is mathematically elegant, but tends to 
obscure the physical meaning of the non-gradient terms 
F n . 
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FIG. 2: Non-monotonic vortex interaction in a system with 
a passive band (i.e. superconductivity is induced in the sec- 
ond band by an intrinsic proximity effect). In limit of zero 
Josephson coupling, the active band would have k = 8k c , and 
thus would be deep into the type-II region. This figure shows 
that a perturbation in the form of a weak Josephson coupling 
to passive band in this case produces a minimum in the in- 
tervortex potential at a very large distance from the vortex 
center. 



IV. NUMERICAL SOLUTION OF THE 
NONLINEAR PROBLEM. 

The linear analysis presented in the previous section 
can only provide information about the asymptotic tail 
of the intervortcx interaction. To determine the actual 
full intervortex potential, especially in the case of strong 
interband coupling, it is necessary to treat the full non- 
linear Ginzburg-Landau theory, something which is not 
possible analytically. In this section we present interac- 
tion energies for vortex pairs (including cases of relatively 
strong interband coupling), computed numerically using 
a local relaxation method. The numerical method which 
we use is the following: A lattice approximant of the 
energy is minimized with respect all the degrees of free- 
dom in the full Ginzburg-Landau functional subject to 
the constraint that vortex positions remain fixed. This 
gives the intervortex interaction energy as a function of 
inter vortex separation. We used high resolutions grids, 
with the number of data points ranging from 1600 x 1600 
to 2400 x 1700 and relaxed each configuration for 50 — 100 
hours on an eight core cluster node. 

In this section the length scale is given in units of 
2//ti, where, as in section III /ti denotes the mass 
of the field | tpi | in the absence of interband couping 
(?/! = 772 = v = 0). Alternatively, the unit of length 
is \/2£i, where £1 = -\/2/Ai is the coherence length of the 
first band in the uncoupled case. Recall that £1 cannot 
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FIG. 3: Intervortex interaction potential for a set of systems 
with two active bands. The systems share the parameters 
given in the table 'common parameters'. The green curve (4) 
corresponds to the case where the bands are coupled by the 
vector potential only. In this case, the ratio of the coherence 
lengths is ^2/Ci = 4. The curve (2) shows the effect of the 
addition of Josephson term, the curve (5) shows the effect 
of addition of mixed gradient term. The curves (1) and (3) 
show the effect of the presence of both mixed gradient and 
Josephson terms with similar and opposite signs. 
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FIG. 4: Gradient, magnetic and potential energy contribu- 
tions to the vortex interaction energy for the parameter set 
corresponding to the curve (3) of the Fig. [3] (black curve). 
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FIG. 5: Non monotonic vortex interaction in systems with 
two active bands and significant disparity in length scales as- 
sociated with the density variation. In the absence of inter 
band coupling (curve 1), the ratio of the coherence length is 
= 12. 
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FIG. 6: Comparison of how Josephson coupling (1) and high 
order density coupling (2-4) affects vortex interaction energy. 
In (3), 772 is chosen so that the densities are approximately 
the same as in (1). In (4), 772 is chosen to give the same con- 
densation energy as (1). Both these parameter values give 
larger vortex binding energy. Similar binding energy as (1) is 
acquired for a significantly smaller 772 (2) . The large condensa- 
tion energy associated with Josephson coupling is responsible 
for the shorter interaction range in (1). 



be identified with physical coherence length when inter- 
band coupling is present. We also measure condensate 
density \tp a \ in units of ui, the vacuum expectation value 
of in the uncoupled case. As shown in appendix [Aj 
this amounts to using scale freedom to set a\ = — 1 and 
pi = 1. In the single component limit, the parameter 
e can then be interpreted as an inverse GL parameter. 
More precisely, n — V2/e, so in these units, in single- 
component limit, the critical value of e is e c = 2. The 
inter vortex interaction energy is given in units of total 
vortex energy, i.e. 2E V where E v is the energy of a single 
isolated vortex. All energies are measured relative to the 
uniform Meissner state (tp a — u a , A — 0, in the notation 
of section III ) . 



A. Weak Josephson coupling to a passive band 

Let us consider a strongly type-II single-component su- 
perconductor, and see how vortex interaction in this sys- 
tem is modified by a weak Josephson coupling to a passive 
band (i.e. a band which has no superconductivity of its 
own which in the context of GL theory manifests itself 
as a positive coefficient 02). 

The Fig. [2] shows the vortex interaction energy in such 
a system. In the limit of decoupled bands, the parameter 
e is here e c /8. Therefore, for zero Josephson coupling 
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FIG. 7: Cross sections of two interacting vortices. The densities are given as 1 — |i/>j| 2 /(|i/i| 2 }, i.e. the condensates are normalized 
to 1 and turned upside down. The magnetic field is also normalized. The left column corresponds to curve 3) of figure ([5| and 
the right column correspond to curve 4) of the same figure. The inter-vortex separation is 3.0 in the upper row, giving repulsion, 
5.4 in the mid row corresponding to the minimum energy, and 7.8 in the bottom row giving attraction. The curves 3 and 4 in 
fig [5] are the ones showing the largest binding energy, a property resulting from them having the largest vacuum expectation 
values in the second band (0.526, 0.5). A careful comparison shows that the recovery of the second band is slower in the right 
column where there is zero inter band coupling. This is a very generic result, Josephson coupling shrinks the disparity in length 
scales. The result of this effect is also seen when comparing the vortex interaction energies of the two system, the curve 3) 
reveals a shorter vortex interaction range than curve 4). 



this system would have k sa 5.7, putting it far into the 
type-II region. Weak Josephson coupling changes the 
length scales as discussed in the previous section and 
adds a qualitatively new feature: the inter vortex po- 
tential acquires a minimum, occurring at a separation of 
approximately 24-\/2£i. 



B. Effects of Josephson and mixed gradient terms 
in case of two active bands 

The figure [3] illustrates the effect of the mixed gradient 
term, as well as of Josephson coupling in a system with 
two active bands. The green curve (4) corresponds to two 
independent bands, interacting only trough the magnetic 
field. The curve (2) corresponds to a system with added 
Josephson coupling which increases the binding energy, 
but decreases the distance where the energy minimum is 
located and slightly reduces the range of interaction. 
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FIG. 8: Transition from type-1.5 to type-II region in a system 
with two active bands. The binding energies are 0.0109 in the 
first curve and 0.0023 in the second curve. The third curve 
give monotonic repulsion. The small binding energy in the 
second case signals that in this particular example the critical 
value for e is close to 0.35. 



The inclusion of a mixed gradient term (shown as curve 
(5)) here has a similar effect on phase difference as the 
Josephson term. When the phases are locked 81—62 = 0, 
effectively this term gives a negative contribution to the 
energy associated with co-directed currents. Thus, for 
this choice of the sign of v, the mixed gradient term also 
prefers phase locking 6\ — 8 2 = 0. 

Due to symmetry, changing rji — > —771 and v — > —v 
does not qualitatively change the behavior of the system, 
as this only results in phase locking with tt difference 
instead. While Josephson coupling increases the energy 
of vortices, and mixed gradients decreases it, their effect 
on interaction energy is the opposite. The decomposition 
of vortex interaction energy into a set of contributions 
from different terms given in Fig [4] illustrates why mixed 
gradients in this case increase repulsion. 

In contrast the curve (1) (blue curve), corresponds to 
the case where v and r\\ have different signs, and so 
there is competition between the gradient mixing and the 
Josephson term with regard to the preferred phase dif- 
ference. The mixed gradient term is minimal for a phase 
locking where 8\ — 62 — n, while the Josephson term is 
minimal for 61 — 9 2 = 0. The result in these simulations 
was that the phase locking was determined by the domi- 
nating Josephson coupling, and that the gradient mixing 
resulted in increased cost for co-directed currents. This 
was the most energetically expensive vortex, but also ex- 
hibited the smallest inter vortex interaction energy. 



C. Solutions with large disparity in the 
characteristic length scales 

Figure [5] shows a set of simulations done with two ac- 
tive bands and a larger disparity in characteristic length 
scales. We start with case when the condensates interact 
only through the magnetic field (blue curve (1)), and 
the density in the second band is 1/4 of the density 
in the first band and the coherence length ratio being 
= 12 (so in the notation of section [TTTj u\ = 4{(2 
and /ti = 12/t 2 ). This allows non-monotonic interac- 
tion to occur at smaller e than above - here we simulate 
at e = 0.7. This gives the smallest binding energy in 
this set of simulations. Adding a Josephson coupling of 
r]i = 0.00556 [shown on curve (2)] gives a substantially 
higher density in the second band, and thus a stronger 
binding energy. Adding a stronger Josephson coupling 
[shown on the curve (3)] gives even larger density in the 
second band. It also shows some of the qualitative differ- 
ences associated with this coupling. Namely, the binding 
occurs at a smaller separation, and the range of the in- 
teraction decreases. It is clearly visible that curve (3) 
crosses the other curves. The Josephson coupling causes 
the second condensates to recover faster (as follows also 
from the linear theory presented in the previous section) , 
and thereby decreasing the range of attractive interac- 
tion. 

Decreasing (3 2 raises the density of condensate in the 
second band. In curve (4), we have reduced it by a factor 
2, thereby increasing the vacuum expectation value of 
the density in the second band by a factor 2. This does 
however not change the length scale, as £ is independent 
of f3, but it does increase the energetic benefits of core 
overlap in the second component. This case shows the 
largest binding energy in this set of simulations. 

Finally, we consider the effect of a higher order density- 
density coupling r]2 between the condensates which is 
shown on the curve (5). The parameter choice here gives 
approximately the same densities as (2) but smaller con- 
densation energy. This should generally make the sys- 
tem (5) recover slower than system(2), resulting in longer 
range of the interaction. A more systematic comparison 
of Josephson coupling and higher order density coupling 
supporting this conclusion is given in Fig. [6] 

Figure [7] displays cross sections of the condensate den- 
sities and magnetic flux in the systems 3-4 in Fig. [5] 
clearly illustrating the mechanism by which type-1.5 su- 
perconductivity appears. 

Let us consider a different example of the appearance 
of the type-1.5 regime in the case where there is a sub- 
stantial disparity in the dominant length scales associ- 
ated with the variations of densities and magnetic field 
penetration length. 

Again, our starting point is a reference case of two 
bands coupled only by vector potential where we choose 
012 so that the disparity in coherence length in absence 
of inter band coupling is £2/^1 = 12. 

Now, we take $2 to be 0.0139 i.e. the same as in curve 
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FIG. 9: Cross sections plots showing condensate density and 
magnetic flux in type-1.5 superconductors with small e. a) 
Curve (1) in Fig. [8] at a separation of 7.2, corresponding to 
the energy minimum, b) The same system at a separation of 
9.6. c) Curve (2) of the same figure at a separation of 9.6, cor- 
responding to the energy minimum. The three distinct length 
scales associated with two band superconductors are indeed 
visible. Despite the inter-band Josephson coupling, there is 
a significant disparity in the recovery of the two condensates 
in all the plots. The third length scale, penetration depth, 
visibly differs between the two systems, and is responsible for 
the differences in inter vortex interaction. 



(4) of Fig. [5] and choose the Josephson coupling to be 
771 = 0.00556. Then, we successively decrease e to see 
where the system crosses over from type-1.5 to type-II. 
The outcome can be seen in Fig. [5J The first curve gives 
a binding energy of 0.011, the second gives 0.0023 and 
the third curve shows system crossing over into type-II 
regime by showing monotonic repulsion. Given the small 
binding energy of the second curve, the cross over from 
type-1.5 to type-II, in this example, is close to e = 0.35. 

The crossection plots of two of these systems given in 
Fig [9] illustrate how the system crosses over from type- 
1.5 to type-II as e is decreased resulting in dominance of 
the repulsive force originating in the electromagnetic and 
current-current interaction. 



V. CONCLUSIONS 

In this paper we presented an analytic and numerical 
study of the appearance of type-1.5 superconductivity in 
the case of two bands with various kinds of substantial 
interband couplings. In all the cases which we considered 
we demonstrated that the system possessed three funda- 
mental length scales: one length scale I /ha associated 
with London magnetic field penetration length while the 
other two fundamental scales l/fJ-1,2 ar e associated with 
characteristic length scales controlling variations of den- 
sity fields. In the limit of two condensates coupled only 
electromagnetically the length scales l/^i,2 are associ- 
ated with independent coherence lengths of two fields. 
However we show that introducing a nonzero Josephson 
and quartic density-density couplings makes both den- 
sity fields decay according to the same exponential law at 
very large distances from the core while, at the same time 
the system still possesses two fundamental length scales 
which are associated now with variation of linear com- 
binations of density fields rotated by a "mixing angle". 
The third fundamental length scale in that regime is the 
London penetration length and thus two-band systems 
with couplings allow a well defined type-1.5 behavior. 
Next we studied the effect of mixed gradient terms and 
showed how the type-1.5 regime is described in that case. 
We showed that in the case of a substantial mixed gra- 
dient coupling the definition of three fundamental length 
scales requires additional care because it produces mode 
mixing which cannot be described by a single mixing an- 
gle. Importantly, we demonstrated that mixed gradient 
coupling can enhance the disparity of the characteristic 
length scales of the density variations. An analogy can be 
drawn between this mechanism and the Seesaw mecha- 
nism in neutrino physics. In the second part of the paper 
we presented a comparative numerical study of type-1.5 
vortices in the different regimes with various intercom- 
ponent couplings. The results were demonstrated in the 
framework of a two-component Ginzburg-Landau model 
with local electrodynamics. However we expect that the 
described type-1.5 behavior is similarly present in lower- 
temperature regimes and in two-component models with 
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non-local electrodynamics. 

The concept of type-1.5 superconductivity can be 
straightforwardly generalized to N-component case. 
There it can occur in systems where characteristic length 
scales are £1, < A < £fc+i, ...,£/v and there are ther- 
modynamically stable vortices with nonmonotonic inter- 
action. 

Besides multi-band superconductors and coexistent 
electronic and nuclear superconductors our model can be 
realized in artificially fabricated layered structures made 
of type-II and type-I materials where one can control and 
tune intercomponent Josephson coupling. 

Note added: After the completion of this work, it was 
verified in a microscopic calculation which does not in- 
volve (1— T/T c ) expansion, that the GL model ^ should 
quite accurately describe the vortex physics in two-band 
superconductors in a quite wide range of parameters and 
temperatures^. 



Let us define the rescaled quantities 
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Appendix A: Units 

In this section we give an explicit mapping from our 
representation of the GL model to a more common text- 
book representation. Consider the Ginzburg-Landau 
model in the following quite usual units 



2mi 

h 2 

2m 2 
-vh 2 Re 

+j l(V x Af 

+ai|?/>i| 2 + - 



{(V-i^A)^-(V + i^A)^*} 



a 2 \ip2\ 2 + - 



-mlV'ill^l cos(6»i 



vM 2 \^\ 2 (Al) 



F = 



(V + ie* : A)ipi 



(V + ie*A)ip 2 



-vRe 
1 



| (V + ie*A)ip 1 ■ (V - ie*A)ip2 } 
V x A\' 2 



+m\^iU2\cos(e 1 -e 2 ) + m\H 2 \H 2 , (A3) 

which, on dropping the tildes, coincides with the repre- 
sentation ([2]) used in this paper. 

Throughout the paper, it is assumed that band 1 is 
active, that is, ct\ < 0. It is convenient to rescale the 
expression ^ for F further so that a\ is normalized to 
— 1 and /3i is normalized to 1. This can be achieved as 
follows. Recall (see section III A | that in the absence of 
interband couplings (i.e. when T]i = 772 = v = 0) conden- 
sate 1 has decay length-scale = (—4a!) -1 / 2 . This 
scale is more usually specified by the coherence length 



£1 = 



V2 



a/^2^ 



(A4) 



We emphasize once more that, in the presence of inter- 
band couplings, £1 in not the coherence length of conden- 
sate 1. This is the purpose of the hat, to remind us that 
this is a genuine coherence length only in the uncoupled 
case. Recall also that the vacuum density of condensate 
1 in the uncoupled model is 



"1 



-ai 



(A5) 



Our second rescaling amounts to using v^i as the unit 
of length and Hi as the unit of condensate density (along 
with compensating rescalings of F, e* and A). Explicitly, 
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Substituting these into ^ yields 



F 



- \(V + ieA)^] 2 + ^ 



(V + ieA)ip 2 



-vRe { ( V + ieA)ipi ■ (V - ieA) ip* 2 } 
-i|Vx A| 2 

-I^l| 2 + ^l^l| 2 -«2|^2| 2 + f |^ 2 | 2 

-jji|^i||^2| cos(0i - 2 ) + ?72|^i| 2 |^2| 2 . (A7) 



This (having dropped the bars) is the GL energy used in 
section |IV| for the purposes of numerical simulation. 

Finally, we note that the single component GL model 
obtained from (A7) by setting ^2 = has penetration 



depth A = l//i a = 1/e and coherence length £ = l/v2, 
and hence GL parameter 



k = A/C = 



V2 



(A8) 



So, in the parameterization used in section |TVj one may 
regard e as an inverse GL parameter for the associated 
single band model. The value of e corresponding to the 
Bogomolny limit of one-component theory is e c = 2 in 
this interpretation. 



Appendix B: The spectrum of H 

Let T-L be a real, symmetric 2x2 matrix both of whose 
eigenvalues are positive, and let H(v) = P(v)^ 1 'H where 
P{v) is defined in equation ( |56| ). We wish to show that 
the eigenvalues \\{v),\2{v) of T-L(v) are also real and 



positive for all < v < 1. First note that H(0) = H, 
so the conclusion holds at v — 0. Further, \ a (y) depend 
continuously on v. Now det'H(^) = det"H/(l — v 2 ) ^ 
0, so neither eigenvalue ever vanishes. So \ a (v) remain 
real and positive, unless, for some v = v* € (0,1), they 
coalesce and bifurcate into a complex conjugate pair. But 
then, at v = we have Ai(i^*) = A2(^*) = A* G (0, 00) 
and hence "H(t^) = A*Ia. But then 

H{u) = PM-^)^,) = A,F(i^)- 1 P(^) (Bl) 

which is symmetric, and hence has real eigenvalues, for 
all v € (0, 1). Hence a bifurcation to a complex conjugate 
pair of eigenvalues is not possible, and we conclude that 
%{v) has real, positive eigenvalues for all v € [0, 1). 



Appendix C: Calculation of long range inter-vortex 
forces from the linear field asymptotics. 

Here we outline how asymptotic intervortex forces are 
calculated from the linearized asymptotic behavior of the 
fields. In the above we use a two-component generaliza- 
tion of the method previously developed by one of us in 
the context of the single component GL model 19 . The 
key idea is to identify the vortices with static topological 
solitons in a relativistic extension of the TCGL model to 
2 + 1 dimensional Minkowski space, which could be called 
a two component abelian Higgs model. Viewed from afar, 
the solitons in this theory are identical to the fields in- 
duced by suitable point sources in the linearization of the 
model, so the forces between well-separated solitons ap- 
proach those between the corresponding fictitious point 
particles, mediated by linear fields. The nature (attrac- 
tive or repulsive) and range of such forces can then be 
computed. In this appendix we take the opportunity to 
explain the method in the simplest possible context, with 
the aim of making it more transparent to a wide reader- 
ship. Despite being simple and pedagogically motivated, 
the calculation we present is, as far as we are aware, new, 
although the result itself has been derived previously by 
other means. 

Consider the sine-Gordon model, which consists of a 
single scalar field <f> in 1 + 1 dimensional Minkowski space, 
evolving according to the Euler-Lagrange equation for 
the action S = J Cdt dx with Lagrangian density 

1 

2 ( 

It is useful to think of <p as an angular variable, with 
period 27r, so that the model has a single ground state 
(or "vacuum"), = 0. It also has static kink (fix solutions 
in which (f> tends to the vacuum as x — > ±00 but winds 
once anti-clockwise. Explicitly, 



C 



(CI) 



4>k(x) — 4 tan 



-1 e x-x 



(C2) 



where xq is a free parameter. These are topological soli- 
tons (topologically stable, spatially localized lumps of en- 
ergy) analogous to the vortices of the GL model. Their 
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energy density £ - 2\m 
a lump centred at x — xq 
asymptotic form 



cos fa is localized in 
At large \x\ the kink has 



4> K (x) 



£ e -\x-x \ 
2tt - 4e-\ x - x °\ 



x — > — oo 
x — > oo. 



(C3) 



We wish to identify this with the field induced by a suit- 
able point source in the linearization of the field theory 
about the ground state <fi = = 2ir. The linearized field 
theory has Lagrangian density 



Co = \d^W4> 



1 



(C4) 



obtained by expanding C about 4> — to quadratic order. 
The corresponding Euler-Lagrange equation is the Klein- 
Gordon equation fa t — <j) xx + = for a scalar field of 
mass 1, whose general static solution is 



cie 



c 2 e 



(C5) 



None of these reproduces the asymptotic kink on the 
whole real line: one needs c 2 = for x < xq and c\ = 
for x > xq. To reproduce the kink's asymptotics we must 
introduce a source term into Cq, 

C ^ C + K(j> (C6) 

where k(x) is the source. The field equation is now 

fat -<l>xx + <f> = k(x). (C7) 

If we take n to be a scalar monopole of charge q placed 
at x = 0, that is, k(x) = q6(x), the induced fie ld is 



-|*| 



We deduce that the asymptotic kink ( C3 ) 



will be induced by the point source 

kk{x) — mS'(x — Xq), m = 



(C8) 



which may be interpreted as a scalar dipole of moment m 
(here 6' denotes the derivative of the delta distribution) . 

So, viewed from afar, the kink soliton is identical to a 
scalar dipole in the linear theory. On physical grounds, 
the interaction energy of a pair of kinks held a fixed dis- 
tance apart should therefore approach that between a 
pair of scalar dipoles as the distance grows large. The 
latter interaction energy is easily computed. Consider 
the field 4> induced by a static source k(x) which is itself 
the sum of two static sources Ki(x) + K2(x). Since the 
field theory is linear, 4> = fa + fa, where fa is the field 
induced by ft, . The total action of the field 4> and source 
k is 



S = 



WW 



+ <h)d^((h + fa)-\(fa + fa) 



+ + K2){4>1 + 4>2)j dx dt 

= S 1 + S 2 + J n^dxdt (C9) 
where Si is the action of (fa, and we have integrated 



by parts and used the fact that fa satisfies (C7| with 



source From this we extract the interaction La- 

grangian for the pair of sources Kx,k 2 , 



Kifa dx. 



(CIO) 



The case of interest is where Ki are scalar dipole sources 
of moment m, located at Xi. One then has k,\(x) — 
mi5'(x — x\) and fa(x) = \ m ii^ e ^ x ~ x ' 2 \ an d so 



J int 



miS'(x — Xi)fa(x) dx 
dip2 
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x—xi 
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Since the interaction Lagrangian depends only on the 
dipoles' positions, we can identify —Li nt as the interac- 
tion potential, 



1 



-mi?7i2e 



-\xi— sj| 



(C12) 



From this we see that like scalar dipoles (mi,rri2 same 
sign) repel one another, while unlike dipoles (mi,m 2 op- 
posite signs) attract. In the case of kinks, mi = m 2 = 8, 
so 



V KK = 32e-l Xl - a:2 l 



(C13) 



in exact agreement with a formula of Perring and Skyrmc, 
which they found using a direct superposition ansatz^Sl. 
Kink-antikink interactions can also be recovered from 
Since the antikink is just a reflected kink, 
4>k(—x), it coincides asymptotically with the 



(C12) 



field induced by a dipole of moment m — —8, so 
Vkk — ~ 32e~' Xl ~ a:2 ' , again in agreement with Perring 
and Skyrme. So well-separated kinks repel one another, 
while kinks and antikinks attract one another. 

The same basic method works for vortices in the TCGL 
model, though the details are more complicated. One 
must linearize the TCGL model about the ground state in 
real ipi gauge, but now there are three fields rather than 
one, A, ipi and ^2, and ipi,ip2 are (generically) directly 
coupled. The coupling is removed by expanding ipx,tp2 
in a basis of normal modes (eigenvectors of the Hessian 
of the potential). The linearized theory is then identi- 
fied with a pair of uncoupled Klein-Gordon models and 
a Proca (massive vector boson) model for A. The point 
source reproducing the asymptotic vortex is a composite 
with scalar monopole charges for the two Klein-Gordon 
fields and a dipole moment for the vector field A. The 
total interaction energy for a pair of point vortices is the 
sum of three terms, two attractive (the scalar monopole 
interactions) and one repulsive (the vector dipole inter- 
action), which can be read off from the quadratic terms 
in the linear theory's Lagrangian density, as described in 
section Mil 
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